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Abstract 

We derive an expression for the mean square displacement of a particle whose motion 
is governed by a uniform, periodic, quantum multi-baker map. The expression is a 
function of both time, t, and Planck's constant, h, and allows a study of both the 
long time, t — > oo, and semi-classical, H — ► 0, limits taken in either order. We 
evaluate the expression using random matrix theory as well as numerically, and 
observe good agreement between both sets of results. 

The long time limit shows that particle transport is generically ballistic, for any 
fixed value of Planck's constant. However, for fixed times, the semi-classical limit 
leads to diffusion. The mean square displacement for non-zero Planck's constant, 
and finite time, exhibits a crossover from diffusive to ballistic motion, with crossover 
time on the order of the inverse of Planck's constant. 

We argue, that these results are generic for a large class of ID quantum random 
walks, similar to the quantum multi-baker, and that a sufficient condition for dif- 
fusion in the semi-classical limit is classically chaotic dynamics in each cell. Some 
connections between our work and the other literature on quantum random walks 
are discussed. These walks are of some interest in the theory of quantum computa- 
tion. 
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1 Introduction 



It is well known that the quantum properties of classically chaotic systems 
differ considerably from their classical versions [1,2,3]. The quantum dynam- 
ics of a classically chaotic system with a finite number of degrees of freedom 
is considerably more regular than its classical counterpart, as indicated by the 
structures developed by Wigner distributions in phase space for simple sys- 
tems, and the fact that almost all definitions of the quantum mechanical ver- 
sion of the Kolmogorov- Sinai rate of entropy production give the value zero for 
finite quantum systems that are classically chaotic with positive Kolmogorov- 
Sinai entropy [4]. Furthermore, it is often difficult to describe the behavior of 
quantum systems for both long times and for small values of Planck's con- 
stant, i.e. in the semi-classical regime. This difficulty is related to the fact 
that the long time limit and the limit of small Planck constant do not com- 
mute, in addition to analytical problems that often make explicit solutions of 
the equations for quantum systems difficult. Therefore it is of some interest 
to study simple models where both the quantum and classical properties are 
accessible to analytic and simple numerical studies, and where the two limits 
can be studied in detail. One example is the baker map. The classical version 
of this map is tractable analytically [5,6,7], and the quantum version is less 
so [8,9], but easily studied using numerical methods. 

It is even more challenging to study transport problems in detail. In recent 
years some attention has been devoted to connecting macroscopic transport 
properties with microscopic chaos, see e.g. [6,7] and other contributions to 
this volume. We find it of interest to study the quantum signatures of these 
relations and look for their appearance in the semi-classical regime. To gain 
intuition about the change in the character of transport properties from quan- 
tum to classical regime a convenient system to study is the uniform, periodic 
multi-baker map. Our interest in this map is stimulated by the fact that the 
classical version is a chaotic system with transport, in this case, diffusion, so 
it provides a model where one can study transport in the context of chaotic 
dynamics. 

Both the classical [10,11] and quantum [12,13] versions of the multi-baker 
map are based upon the baker map, which in its classical version is an area 
preserving, expanding and contracting transformation of the 2-torus onto it- 
self. The multi-baker map is obtained by considering a two dimensional strip 
of unit height in the y-direction and a segment of the real x-axis, either of 
infinite length in both directions or of finite length with specified boundary 
conditions [10,11]. The multi-baker transformation is a combination of the 
baker map with a translation of points in each unit interval to correspond- 
ing points in the nearest intervals to the right or left according to a well 
defined prescription, to be given below. The quantum versions of these maps 
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are obtained by means of a simple quantization where the x-axis is taken to 
correspond to position space, and the y-direction is taken to correspond to mo- 
mentum space [12,13]. The quantum mechanics is obtained by requiring that 
there be an integral number of quantum states in a unit interval, and that the 
time development of these states correspond to the expanding and contracting 
properties of the classical map. Balazs and Voros [8], and Saraceno [9] were 
able to show that these requirements can be satisfied by means of a unitary 
operation on the quantum states, provided Planck's constant is taken to be the 
inverse of the number of quantum states in the unit square. Our work has been 
devoted to extending this quantum baker map to a quantum multi-baker map 
and examining the transport properties of this quantum map [12,13]. Other 
quantizations of multibaker maps as well as the Kapral-Elskens coupled baker 
model have been studied in [14] 

In this paper we focus on the calculation of the mean square displacement 
(MSD) of a particle whose dynamics is described by a quantum multi-baker 
map. This quantity is of importance for the description of the average mo- 
tion of a particle, since the time dependence of this quantity can distinguish 
between sub-diffusive, diffusive, and super-diffusive motion. A particular ex- 
ample of super-diffusive motion is, of course, ballistic motion characteristic of 
the motion of a free particle. 

In the classical multi-baker, the MSD grows linearly in time, ((Ar) 2 ) = t, char- 
acteristic of diffusion, for all times greater than some microscopic time. The 
quantum version shows an asymptotic ballistic, t 2 , growth for what we will 
call uniform maps, which are translationally invariant. This is essentially the 
same as one finds in models of electron transport in periodic solids, such as the 
Kronig- Penney model [15]. One can find examples of non-translationally in- 
variant multi-bakers where the particle is localized with no asymptotic growth 
of the MSD with time. Here we consider only translationally invariant mod- 
els, and leave the large class of non-translationally invariant models for future 
work. 

The classical multi-baker map [10,11] is a simple model for deterministic dif- 
fusion along a one-dimensional lattice. It can easily be constructed to be iso- 
morphic to a random walk with any given probabilities for jumping to the 
right and to the left along the lattice. The quantum multi-baker map [12,13] 
is also an example of a random walk, but in this case, the walk is quantum 
mechanical and has very different properties from the classical map due to the 
interference of probability amplitudes for possible paths. The subject of quan- 
tum random walks has developed a literature over the past few years. Here we 
will also discuss the connection of the quantum multi-baker with other mod- 
els of quantum walks. We mention now that most of the previous studies by 
other authors have been devoted to the so-called Hadamard walk [16], which 
is a special case of the quantum multi-baker for the largest allowed value of 
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Planck's constant. Thus the multi-baker map studied here represents a gener- 
alization of the Hadamard walk to all allowed values of Planck's constant. 

The plan of this paper is as follows. In Section 2 we present the basic equations 
defining the quantum multi-baker map and compare it to the classical version. 
In Section 3 we discuss the MSD for the quantum map, and in Section 4 show 
how random matrix theory may be used to evaluate the MSD for small values 
of Planck's constant, leaving some technical details for a longer paper [17]. 
In Section 5 we compare the results of random matrix theory with numerical 
studies of the map, which shows that random matrix theory is quite effective 
in reproducing the numerical results and in providing an expression showing 
the crossover from diffusive to ballistic motion for fixed, non-zero values of 
Planck's constant as the time becomes large. In Section 6 we discuss the 
connection between the multi-baker map and other models of quantum random 
walks. We conclude with a summary of results, and a discussion of the possible 
implications and extensions of this work. 



2 The quantum multi-baker map 



We begin with the classical multi-baker map [5,10,11]. It is a two-dimensional 
lattice system where the phase space at each lattice site is a square and the 
dynamics is a combination of transport of the phase space densities to neigh- 
boring cells, which models the free flight, followed by a local baker map evolu- 
tion within a square, which models a collision with a fixed scatterer. That is, 
the multi-baker map, M, is a composition of two maps M = BoT: transport, 
T, of phase points to neighboring cells, given by 



and the baker map, B, which acts on the x, y coordinates of each cell, n, 
separately, according to 



T(n,x,y) = < 



(n+l,x,y), for0<a;<l/2, 
(n — 1, x, y), for 1/2 < x < 1, 



B(n,x,y) 



= < 



(n,2x,y/2), for0<:r<l/2, 
(n,2x- l,(l + y)/2), for 1/2 < x < 1. 
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The combination of these two maps is the multi-baker map which is a time- 
reversible, measure preserving, chaotic transformation, with evolution law 

f (n + l,2x,y/2), for0<:r<l/2, 
[ (n - 1, 2x - 1, (1 + y)/2), for 1/2 < x < 1. 

As mentioned above, this classical map represents a simple area-preserving 
model of simple random walk. It can also be considered a simplified Bernoulli 
map for the motion of light weakly interacting particles in a gas of heavy 
scatterers on regular lattice (periodic Lorentz gas). 

The quantization of the multi-baker map is based upon the known quantiza- 
tion of the baker map on the unit square, carried out by Balazs and Voros, and 
by Saraceno [8,9]. We form the quantum multi-baker by transporting some of 
the quantum states to the next cell on the right and the others to the next cell 
on the left, based upon the way the transformed points are moved to neigh- 
boring cells in the classical multi-baker, and then model the scattering by the 
quantum baker map B acting in every single cell. 

To quantize the baker map, we regard the horizontal direction of the torus 
[0, l) 2 as the position axis, while vertical axis corresponds to the momentum di- 
rection. To obtain the Hilbert space we take the subspace of the wave functions 
on a line whose probability densities, |\l/(a;)| 2 , 2 are periodic in both posi- 
tion and momentum representations, respectively: ty(x+l) = exp(i2ir(p q )^(x), 
ty(p + 1) = exp^Tnpp)^^), where (fi q ,(p p G [0,1) are phases parameterizing 
quantization. The quantization of the baker map requires the phase space 
volume to be an integer multiple of the quantum of action [8,9,18,19,20,21]. 
Therefore the effective Planck constant is h = 1/N, for a baker map on a unit 
torus, where N, an integer, is the dimension of the Hilbert space. The space 
and momentum representations are connected by a discrete Fourier transform 
(Pk\qi) = [GN(<Pq, Vp)]fc,J := N' 1 / 2 exp(— i2nNpkqi). The discrete positions and 
momenta are qi — (I + tp q )/N, Pk = (k + ip p )/N, respectively. We take N 
to be an even integer and label the quantum states in coordinate space by 
j, where j = 0, 1, . . . , TV — 1. The states j = 0, . . . N/2 — 1 are called the 
"left" states, with collective wave function, ty L , while the remaining states, 
j — N/2, . . . , N — 1 are called the "right" states, with wave function, ty R . One 
may think of ^r,l as N/2 component vectors. 

Having constructed the Hilbert space, one looks for a family of unitary prop- 
agators parameterized by iV = 1/h which approach the classical map in semi- 
classical limit. The quantum baker map [8,9,19] is then given by a combination 
of two operations. The first operation takes the N/2 "left" states into momen- 
tum states labeled k = 0, 1, ... , N/2 — 1, called "bottom" momentum states, 
by means of the N/2xN/2 Fourier transform matrix G N / 2 ((p q: ip p ). The "right" 
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coordinate states are transformed to "top" momentum states in the same way. 
Now one has transformed N spatial states into N momentum states in a way 
that mimics the classical baker's map. The final step is to express the N new 
states in coordinate representation by means of the N xN matrix, G^ 1 , which 
takes iV states in the momentum representation to their coordinate represen- 
tation. The full baker transformation, B, on the torus is then given by the 
unitary transformation 



B := GV 1 



G N / 2 



G 



N/2 



(1) 



for even N. Other examples and discussions of issues concerning the quantiza- 
tion of area-preserving maps can be found for instance in [8,9,18,19,20,21,22]. 

Using this transformation as a template we can easily express the quantum 
multi-baker map as a transformation of right and left quantum states in unit 
squares labeled by n ± 1 to right and left states in the unit square labeled by 
n, ^R,L{n). That is, in the position representation the evolution of the wave 
function in the quantum multi-baker is given by the equations 3 

* L (n,t+l) 
y R (n,t + l) 

where G^{n) is the discrete Fourier transform 

(G N (n)) kj = -L e - 2 ™( k +<Pp(n))U+? q (n))/N _ (3) 



G N \n) 



G N/2 {n) 
G N/2 (n) 



* L (n-l,t) 
V R (n + l,t) 



(2) 



In principle, the Fourier transformation matrix Gjy(n) can depend upon the 
cell index n. Here we will consider the case where this matrix is independent 
of the cell index, so that the system of equations is translationally invariant 
from cell to cell. We restrict our attention here to finite systems of length L 
with periodic boundary conditions. This condition induces Bloch states and 
an eventual ballistic motion of particles through the chain. We associate L- 
dimensional Hilbert space with the lattice, and since the internal space is 
C^, the Hilbert space of the system is the tensor product C L <g> C N . We will 

3 Note that the construction we use here differs slightly from the one we studied in 
our previous work [12]. We change it to make the construction consistent with the 
general picture that we now have as well as to make model clearer. The properties 
of the two versions of random model should be the same, and the periodic model is 
unaffected. 
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work with the basis defined by \n, ±,i), where n is the lattice site, ± denotes 
left/right half of the square, and % — 0, . . . , N/2 — 1 denotes the state in given 
half in the position representation. Since the states associated with the left half 
of the cell go one step to the right we write ^f L = accordingly ty R = 
Then the general wave function for the multi-baker chain can be written as 

l*) = EE ( 4 ) 

n e=± 



where |\P e (ra)) := P e (n)\^f), P £ (n) := Y,i\n,e,i)(n,e,i\ is orthogonal projec- 
tion of the state vector onto the "left" or "right" subspace at site n. Thus the 
inner product takes form 

EE <*«(")!*«(«)>■ (s) 

n e=± 



We may think of \^ £ (n)) and (\& e (ra)| as corresponding to iV/2-component 
column and row vectors, respectively, while corresponds to an NL-com- 
ponent vector. Since we work mostly in position basis, we use N/2 dimensional 
vectors ^(ra), whose i-ih components are ^±(n,i) = (n, ±, where the 
index n denotes the lattice site of the cell under consideration, as above, the 
notation ± denotes the particular half of the unit square under consideration, 
the left half indicated by + moves to the cell to the right, and the right half, 
denoted by — , moves to the cell to the left. The index i denotes a particular 
quantum state in the right or left half of the unit cell, specified by ± and 
i = 0, 1, . . . , iV/2 — 1. The action of the quantum multi-baker map on the total 
wave function will be written as the action of an operator M defined by Eq. 
(2), as 

= M|tf). (6) 



The operator M is the Floquet operator for the multi-baker map and will be 
used to determine the time dependence of various observables for this system, 
as 

n(t) = Mt'ftM', (7) 



where Q is any observable of the system. 

The complete specification of the model requires values for the phases ip q , <p p . 
There is a considerable amount of freedom for choosing these phases. As men- 
tioned above, one may choose them to vary from cell to cell. A random vari- 
ation of phases from cell to cell produces a disordered multi-baker map, while 
requiring that the phases have constant values throughout the lattice produces 
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a regular, or uniform multi-baker map which is translationally invariant from 
one cell to the next. Many other choices are possible and can be of interest. 
Here we consider only the translationally invariant case. In numerical calcu- 
lations we use the values for the phases chosen by Balazs and Voros [8] with 
<Pq — <f P — 0; those used by Saraceno [9] ip q = ip p = 1/2, which lead to survival 
of additional classical symmetry in the model, as well as more generic values. 
We will see below that the Saraceno phases lead to a non-generic behavior for 
the MSD. 

Let us consider now the structure of eigenstates of the unitary operator M. 
Since we consider periodic boundary conditions every eigenstate corresponding 
to the eigenvalue A = e lK has a Bloch form 



*(n, ±) = exp(inn)^ ± /VL 



(8) 



where 
at or 



is the normalized eigenstate of a modified quantum baker oper- 



G N/2 e iK 



(9) 



Periodic boundary conditions imply e %K,L = 1, thus we have = 2irk/L, k = 
0, 1, . . . , L — 1. For every we have N eigenstates. We will enumerate the 
eigenstates of M by k, n, with n = 0, . . . , N — 1. Thus is the eigenstate 
given by the ansatz (8) with = 27fk/L, corresponding to the eigenvalue 
e llfn , where phases tp k n are counted for a given k from (including) to 2n, i.e. 

< V?n < < 27T. 



3 The mean square displacement (MSD) in the uniform multi- 
baker map 



In order to formulate our calculation of the MSD for the uniform multi-baker 
map in the most convenient way, we introduce two operators, r, v, which rep- 
resent a coarse position operator and a coarse velocity operator, respectively. 
These operators are coarse in the sense that r simply gives the lattice site 
associated with a particular quantum state, and the coarse grained velocity v 
is given by v = MVM — r. We argue elsewhere [17], that for a translationally 
invariant system v has a very simple form with values ±1, given by the change 
in cell index for each quantum state. Thus 
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r\n, ±, i) = n\n, ±, i) 

v\n, ±, i) — ±|n, ±, i), (10) 
so that explicit expressions for the operators r, v are 

r = n|n, e, i){n, i, e\, 

n,e,i 

v=y~] e\n,e,i)(n,e,i\. (11) 

n,e,i 

Suppose we prepare the system in a pure state |\&). Then the mean square 
displacement of the particle starting in this state is given by 



<(Ar) 2 (t)>* = ((M\M* - r)% = ((£ ^) 2 >* = £ (^>*, (12) 

r=0 Tl,T2=0 

where (A)* := = Tr (|^)(^|A). Depending on the original state we 

have a distribution of possible results. To characterize it we can calculate its 
average over all the possible initial states, and the root mean square deviation 
from the average, which quantifies the spread of the results, or the quality of 
the prediction based on the average. In this Section we find the expressions 
for both the equilibrium MSD as well as for the equilibrium fluctuations of 
this function. In the next section we approximate the average results using 
random matrix theory, and then we compare them with numerical evaluation 
of the exact formulas. 

Since the average over all the possible pure states gives the most incoherent 
mixture, we obtain a simple expression for the equilibrium MSD as 



<(Ar) 2 (t)>= E (v n v n ), (13) 
n,T2=0 



where 



(A) := Tr (g cq A) = -^Tr (A), (14) 



L is the length of the chain (assume periodic boundary conditions), and iV is 
the dimension of the Hilbert space for a single cell. Clearly, £> eq = 1nl/(LN) 
is the equilibrium state for all the QMB, that is, it is time invariant 

M&qM 1 = g cq (15) 
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and it maximizes the entropy S = — Tr (glng). Time invariance implies that 
the velocity autocorrelation function 

C T1 , T2 := (v Tl v T2 ) = (M tTl t;M Tl M tT2 t;M T2 ) (16) 
satisfies C TliT2 = C n _ T2) o = Co )T1 _ T2 = C Tl _ T2 . Thus we can write 

((Ar) 2 (t)) = ]TQ, + 2 £ C 71 _ 7a = C t + 2 2(t-r)C T . (17) 

r=0 Tl>T2=0 T=l 

This leads to a Green-Kubo like formula for the diffusion coefficient 

1 -I OO OO 1 

D = lim -((Arm) = ? C + £ C T = £ ^ - A 

Zt 2 r=l r=0 2 

paralleling the classical case. However, here the diffusion coefficient turns out 
to be infinite. 

Thus the MSD is written as the sum of time correlations of the velocity, just 
as in the classical case, but the difference in mechanics will lead to important 
differences in the time development of the MSD. 

The translational invariance of the system, together with the uniform equilib- 
rium density, imply that instead of summing over all the N quantum states in 
the L cells, we can sum only over the N quantum states in one cell, see Ref. 
[17]. Thus we can write 

LN 1 J 

N i,e=± 

= T? ^ e ( n = °' ^l MtT ^ MT | n = 0, e, i). (18) 

N i,e=± 

Thus we have to sum, with appropriate signs, over all the possible cyclic paths 
starting in a fixed cell. The next simplification provided by the translational 
invariance of the model is that we can replace the unitary operator M for the 
full multi-baker in Eq. (18) by the operator B acting within a unit cell. This 
replacement follows from the observation that the velocity operator is diagonal 
in the representation based on the states in the cells, with a block structure 
that is identical from one cell to the next 4 . Therefore, the velocity operator at 

4 This would not be the case for the disordered multi-baker maps where the quan- 
tization phases vary from one cell to the next. 
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time step t takes on the same value in all of quantum states that are periodic 
images of any state in a fixed cell. Thus we can express the velocity correlation 
function, C T , as 



C T = -^Tr [B tr JB T J] 

N 1 J 



(19) 



where the matrix J is the velocity operator restricted to a single cell, in position 
representation given by 



1jv/2 

— ljy/2 



(20) 



Here 1jv/2 is the N/2 x iV/2 unit operator. 

To proceed further we must now use specific properties of the quantum baker 
map B, and in particular, its spectral properties. Since B is a unitary operator, 
its eigenvalues lie on the unit circle and define a set of N phases, ipj or quasi- 
energies. We will denote the corresponding eigenstates by the Dirac kets 
Thus the spectral problem for the operator B is 



B|j) 



\J)- 



(21) 



Then we can express the velocity correlation function C T as 



C T = -J2e-^(k\3\j)e^(j\3\k} 



— Y |J 1 -fc| 2 e i(w "*' fc)T 



j>k 



The matrix elements of J satisfy 

£ |Jjfc| 2 = Tr J 2 = N = J2 |Jii| 2 + E 21 Ji*| 



(22) 



(23) 



Therefore 



N + 12 \3 jk \ 2 {e iiVj ~ Vk)r + e- i{fP '- VH)T - 2) 
j>k 
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^-I^WM^ (24) 

j>k Z 

We can now substitute these results into Eq. (17) for the mean square dis- 
placement to obtain 

<(Ar) 2 (t)>=t + E(t-r)CV 

T = l 

= t + t(t - 1)1 £ | J,,f + 1 E | J jfc | 2 £ rKe^^)^)(25) 

j j>fc T=l 

The last expression is particularly useful for averaging and will be used in the 
next section. 

The sums over the intermediate times r in Eq. (25) can easily be carried out. 
We assume now that there is no degeneracy in the spectrum of quasi-energies 
and, after some algebra, we obtain an expression for the MSD as 



1 o „; n 2 (<Pj-<Pk)t 

«*fM> =4 T. fel 2 + 1 E fel 2 g -^ (26) 

iV j iV j>fc Sm 2 

j s j n 2 (vj-yfc)* 

^ sm 2^ 

Whenever two eigenphases (pj,(pk are equal, the contribution to the sum 

gin 2 (VjZZkH 

— y.-y. must be replaced by t 2 . 

sin 2 J_ 

We see that there is typically a ballistic contribution coming from diagonal 
and possibly some degenerate terms. The other contributions are oscillatory 
and usually negligible in the long time limit. It can happen, however, that 
the ballistic contribution disappears altogether and we have only oscillations 
which means that the particle is localized. An example is provided in Section 5 
and an explanation in more general context is given in the final Section. It is 
interesting, that the obtained result is completely independent of the length 
of the system. This is not true for the fluctuations of the MSD. 



It is not easy to evaluate this expression, Eq. (27), analytically. Therefore in 
section 4 we present an approximation based on the random matrix theory, 
and in Section 5 we evaluate it numerically for some particular choices of 
parameters, and compare it with the RMT averages. 
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3.1 Fluctuations in the MSD 



We have computed the MSD as an equilibrium average where all eigenstates 
have the same weight. One might ask about the dependence of the average 
square displacement for an individual quantum state. In order to say some- 
thing about the fluctuations of the average square displacement from one state 
to the next, we now consider, A*r 2 (t) the average square displacement for the 
quantum state |\l/), where A*r 2 (t) = Ar 2 (t)\^f). We will try to characterize 
the fluctuations of the square displacement among various quantum states by 
calculating, in so far as it is possible, the mean square fluctuation of A*r 2 (t), 
which we denote by A 2 (t), where 

A 2 (t) = /[A*r 2 (t)-(Ar 2 (t))l 2 \ , (28) 



The average (A(^f))^ is defined as follows: if |\P a ) is an arbitrary orthonormal 
basis and b a := (^ a |^) are the complex coefficients of expansion of |^) in this 
basis then 

/ A(\u\\ - ^ ■ ■■d 2 b LN 5{l - El L ! \b a \ 2 )A&) JcM 2LN A(y) 

K^yvj/y-— ,, 9 , jo, xn ^nl il .ox - rjn • y zy > 



Thus it is an average over 2LiV-dimensional sphere. 

Before we calculate the right hand side of Eq. (28), we note that the equilib- 
rium mean square displacement can be written in a very similar way as an 
average over the space of all the states 

<(Ar(t)) 2 (t)> = <A*r 2 (t)>*. (30) 



Indeed, using Eq. (29) we obtain 

<(Ar(t)) 2 (t)> = E<W(^|A^)I^> (31) 

a,/3 



Only contributions from a — (3 survive due to phase averaging, and the diag- 
onal terms carry the same contribution, therefore (6*6/3)* = 5 a ,f3l \NL). If we 
insert this relation into the right hand side of Eq. (31), we recover Eq. (14) 
with the quantity A = Ar 2 (t). 

Returning to Eq. (28), we see that the right hand side can be written as 
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t-1 LN 

A2 W= E E (KbK b M^\ V ^ V rn\^p){%\VkVl\^ S ) 

n,m,k,l=0 a, fi, 7, 5=1 

-[(Ar 2 (i))] 2 . (32) 

Observe that if at least one phase, or equivalently one subscript, is unpaired, 
the average (6*6/36*65) will be zero. Therefore the average of the four coeffi- 
cients is non-zero only if the indices are paired. This can be arranged in three 
ways, so that 

(6*6,36*6,5) = Ci<5 a/3 <5 75 (1 - 5 fa) + ci5 a5 5 7/3 (l - 8 a p) + c 2 <!>a/A<5 ( W, (33) 
where 

Ci = (N 2 N 2 > (34) 
c 2 = (|6 1 | 4 ) (35) 

Let A = Em^oVmVn- Then 

A 2 (t)= Cl £<*a|4*a><*0|4*0> + Cl 2<*a|^|^>(*/3|A|* a ) 
+ ^C 2 (* a |A|* Q )(* Q |A|* Q )- ' 



1 -T 4 



l^a>| 2 - 


[ 1 Tr^l 







(.36) 



Since the term J2a \(&a\A\ty a )\ 2 is basis dependent we must have c 2 = 2c±. 
Also, 

EI6«l 2 N 2 = (EI6«l 2 N 2 > = i- (37) 

Since there are LA" quantum states, this reduces to 

c 2 LN + Cl LN(LN - 1) = 1, (38) 

therefore 

1 2 
Cl ~ LN{LN + 1) ° 2 ~ LN(LN + 1)' (39) 

To proceed further we need to use the results obtained earlier following Eq. 
(9) for the eigenstates and eigenphases of the multi-baker operator. The eigen- 
states have the form of periodic Bloch states and the eigenphases appear in L 
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sets of N eigenphases. After some lengthy, but straightforward, calculations 
one finds the final result for the fluctuations of the average square deviations 
as [17] 



A 2 (t) = ci[(Tr A) 2 + Tr A 2 



L 2 N'- 



■(Tr Af 




sin sin 



sin p 2 



sin 



(40) 



The main source for the system size, L, dependence of this fluctuation result is 
the factor of L^ 1 in front of the whole expression. This provides the usual L~ l l 2 
decay of the root mean square fluctuation which suggests that in the large 
systems the equilibrium average is typical for almost every initial, pure state. 
The eigenphases y?" also depend on the system size, so that in the absence 
of a further evaluation of the sums, one cannot precisely determine the size 
dependence of A 2 (t), however, it is reasonable to expect that the contributions 
for different Bloch vectors are of the same order. Numerical results suggest that 
this expectation is correct for the quantum multi-bakers, and the final result 
for the right hand side of Eq. (40) decays inversely with L. This is shown in 
Figure 1. The left picture shows the time and size dependence of the relative 

fluctuations j&^jj on a log- log plot. The right picture shows the same data 
with each curve multiplied by square root of L. The scaling is remarkable. 
Plots were obtained through numerical evaluation of exact formula (40) for 
N = 50, ipg = ipp = 0. Using the methods of random matrix theory discussed 
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Fig. 1. Time and length dependence of the relative fluctuations ^r^(t)) • ^e left 
picture shows log-log plot of the unsealed data, the right picture shows each curve 
multiplied by square root of L. Plots made for N = 50, ip q = ip p = 0. 



in the next section, in the limit of large N, one can replace the average over 
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the Bloch states, labeled by the subscript a, by averages over random matrix 
ensembles. We have not yet done this and leave it for future work. 



4 Random matrix theory for the mean square displacement 

To proceed further we need to know the quasi-energies and eigenvectors of 
the operator B. For large N it is not possible to determine these quantities 
by analytical means. Instead one uses numerical methods to diagonalize B. 
The results obtained this way will be discussed in the next section. Here 
we will show how random matrix theory [2,23] may be used to evaluate the 
approximate value of the MSD. In a sense, the use of random matrix theory 
must be considered to be a kind of uncontrolled approximation method since 
there are no theorems indicating that the distributions of quasi-energies are 
described by any of the three random matrix ensembles, and there are no 
similar results for the eigenvectors. However, since the classical baker map 
is an example of a strongly chaotic system, we can expect that one of the 
two ensembles, the circular orthogonal ensemble, COE, or the circular unitary 
ensemble, CUE, might provide a very good approximation when it is used to 
evaluate the right hand sides of Eq. (27), see ref. [9,24]. The results obtained 
this way will, of course, be averages over the ensemble of random matrices, and 
not characteristic of any one matrix in the ensemble, and perhaps not of the 
baker map itself. We will see shortly that the numerical results for the baker 
map are in very good agreement with the results of random matrix theory, 
so the connection between the use of random matrix methods for classically 
chaotic systems is supported by this work. 

We begin the application of random matrix theory by determining the average 
values of the moduli of the matrix elements Jjk — ( j | J ] k ) where \j), \ k) 
indicate eigenvectors of the unitary baker operator B. We begin with the 
diagonal elements, which we first express in the position representation of the 
N quantum states defined in the baker cell. That is, we write 

/JV/2-1 JV-1 \ 2 

\J»\ 2 = MJ\j)\ 2 = [ E K«b)l 2 - E l(«b')l 2 ' ( 41 ) 

\ a=0 a=N/2 ) 

where the kets \a) are basis vectors in position space. Since the eigenkets are 

normalized, that is, J2a=o \( a \j)\ 2 = 1> we ma J se ^ s j = Ea=o 1 l( a |j)| 2 to 
obtain 

|J ii | 2 = (2a i -l) 2 = 4 a ;-4a J - + l. (42) 
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To get the random matrix theory average value, we observe that the joint 
probability of single eigenstate components in the CUE and COE ensembles 
are given by [2,25] 



2M-1 



P({*}) = 05(1 - |*| 2 ) = c5(l - J2 W 2 )- (43) 



a=0 



In the COE ensemble the eigenstates can be chosen real which implies that 
P({^}) is a uniform distribution on iV-dimensional unit sphere, and M = 
N/2. In CUE ensemble the eigenstates are complex implying that P({^/}) is 
uniform on a 2iV-dimensional unit sphere {\^j\ 2 = {^j) 2 + (O^) 2 ), and 
M = N. These averages can be calculated in straightforward ways and one 
obtains (s) = 1/2 for both ensembles. However the averages (s 2 ) differ in the 
two ensembles. We obtain 

= r(f + 2) r(M) = m + 2 
^ I r(f) r(M + 2) 4(m + i)' 1 ; 

where M = N/k with k = 2 for the COE and k = 1 for the CUE. This leads 
to an evaluation of the averages needed for the diagonal terms, with the result 
that 

(fef ) - J^~ k - (45) 



To get the average of the off-diagonal terms, we use (23) which, after averaging, 
leads to 

(l ^ = (N-ll N + kY (46) 



with k given above. 

4.1 The MSD in the CUE ensemble 

Next we calculate the average of e l( -' fil ~' P2 ^ T in CUE ensemble. To do this we 
need an expression for the pair correlation function, R({pj,{pk) of two angles 
in this ensemble. This is calculated in some detail by Mehta [23] , and we use 
the expression given there. 



2tt 2tt . . 

T R{(p j ,(p k ) 



N(N-l) 
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2tt 2tt 



A 



o o 



4vr 2 (A- 1) 



1 - 



sin 



2 N(ipj-ip k ) 



A 2 sin 



2 (<fij~Vk) 



(47) 



At t = the first term is A/ (A — 1), and it vanishes for r > 0. We calculate 
the second term by converting the angular integral to one in the complex 
plane. By changing the variables to u — ipi — (fi2,v = <f2, and then setting 
z = e m , dz = iz du, we obtain, for r > 0, 



-1 



2tt 



-J 



2vrA(A - 1) 



27riN(N - 1) 



sin 2 ^ 

due™ 2 



sin 



2 u 



dzz T - N g(z), 



(48) 
(49) 



where 



iV-l N-l 2N-2 

g(z) = J2a p z? = (J2 z k f = J2(k + l)z k + J2 (2N-k-l)z h 

p k=0 k=0 k=N 



(50) 



Since g(z) is analytic, then (exp[i(</?i — <p?)t\) = if r — A > 0. Otherwise, 
we can write 



^/^- gW = i 9 «(0), 



(51) 



where p — N — r — 1, and is the p-th derivative of g, and obtain 



i(fi— <P2)t 



-a, 



N(N-iy 



(52) 



Thus, including the contribution from the first term at r = 0, we obtain 



(eMn-v^cuE = < 



1 for r = 

for < r < A, 
for r > A. 



(53) 



We now substitute the results, Eqs. (45), (46), and (53) into the expression 
Eq. (25) the for mean square displacement. This leads to 



t-i 



{(Arf)=t + t{t - l){\Jjj\ 2 ) + (A - l)(\J^ k \ 2 ) £(t - r)(e- + e~ ic 



r=l 
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M-1 



t + 



N + l 



'('-^ V + 1)^-1) 5"-^-^ (54) 



where M = t for t < N, and M = N for t > N. Carrying out the required 
sums we obtain the final result for the MSD in the CUE ensemble: 



<(Ar) 2 (t)> = 



t + 



t{t ~ 1) 
N + l 



t-2 
3(N-1) 



{N + l 



1 2 N 
t + — 



for t < N, 
for t > N. 



(55) 



4.2 The MSD in the COE ensemble 



The calculation of the MSD in the circular orthogonal ensemble proceeds in 
very much the same way as in the CUE, the only difference being in the pair 
correlation function for the quasi-energies, ipj. For the COE the two-point 
correlation function is [see Mehta[23], Eq.(6.2.4), Eq.(10.3.41)] 



i2 2 (tf, ^)=det 



a N (0) a N (-d - if) 
ctjv(v?-i?) <7jv(0) 

= [(^(0)) 2 ] (0) -[(^(^-v.)) 2 ] (0) . 

Here, is a quaternion, and det is the quaternion determinant, with 

[(a N (0 - y?)) 2 ] (0) = (S N (0 - y?)) 2 + DS N (0 - <p) JS N (d - ip). (56) 



These quantities are given by Mehta as 



2tt sin | 2ir k {i 



1 N ~ 1 d 
DSn{$) = ^ E iPk ^ = ^ S n(#) 

1 ^ 1 



eft? 



-E 



(57) 



where p fc = |-f +£;, fc = 0, 1, . . . , JV-1, and q = ±\(N+1), ±|(iV+3), . . .. 



Then 



coe = N ( N _ 1 ) [ S i(t) + S 2 (r) + S 3 (t)}, 



(58) 
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where 



5i(r)= I d& d<pe 



AT2 



- - -i N-l 

S 2 (t) = - J d$ J dipe^-^— J2 e l{pk+Pl) ^\ 



k,l=0 



S 3 (t)= I dd I ^e^-^-L^^e^X*-*'). 



(59) 



p,q 



The evaluation of these sums and integrals is relatively straightforward, and 
leads to 



/ [iar]\ _ ) 

\e L '/c'oe — \ 



1 for r = 

j^l-N + 2r[/(f + r) - /(f)]] for < r < N, (60) 
j^i-N + 2r[/(f + r) - /(r - f )]] for r > TV. 



Here /(T) is defined by 



f(T) := £ 



1 



k=i 



2k -I 



1 1 



(61) 



This function has at most a logarithmic dependence on its upper limit for 
large T. For the MSD in the COE we then obtain 



<(Ar) 2 (t)> 



t + 



t(t - 1) 
N + 2 



1 + 



N 



t-2 
3(N-1) 
N 



[N + 2 3 3(iV + 2) 



+ 5< for t < N, 
+ 5> for t > N. 



(62) 



Here 5 <t> are small corrections to the explicit formulae that have to be evalu- 
ated numerically. Note that the explicit results given here for both ensembles 
can be combined into the single expression 



<(Ar) 2 (t)> 



t + 



t(t - 1) 
N + k 



k-l + 



t-2 
3{N- 1) 



+ (k — l)5 < for t < N, 



k 

. N + k 



r + — - 



N N(k-1) 



(63) 



3(N + k) 



+ {k- 1)5> 



for t > N. 



Those results are shown in Figure 2 for iV = 200. The COE results are two 
close curves, where the higher is the result given in Eq. (63) for k = 2, while 
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in the lower curve the corrections 5 <t > have been neglected. Three asymptotic 
estimates t,t 2 /N,2t 2 /N are also plotted. Inset shows the region t = 100 to 
t = 300 where the differences between the two COE results are most pro- 
nounced. To our delight, RMT average lends classical diffusion as the short- 



1 00000 




time t 

Fig. 2. Log- log plot of the ensemble averages of the mean square displacement using 
RMT. Both CUE and COE results are plotted. The COE results are the two close 
curves, where the higher is the result given in Eq. (63) for A; = 2, while in the 
lower curve the corrections 5<,> have been neglected. Three asymptotic estimates 
t, t 2 /N, 2t 2 /N are also plotted. Inset shows the region t = 100 to t = 300 where the 
differences between the two COE results are most pronounced. 

time prediction in both cases, the CUE average being "more classical". It is 
worth emphasizing that the classical behavior persists up to the Heisenberg 
times ~ h^ 1 = N rather than the Ehrenfest time ~ \rih~ 1 = InN. On the 
other hand, for times longer than the Heisenberg time we observe ballistic mo- 
tion. The ballistic coefficient is proportional to the effective Planck constant 
therefore it disappears in semi-classical limit. 

In other words, fixing the time and performing semi-classical limit (h — > = 
N — > oo) we obtain ((Ar) 2 (t)) = t, which is the classical result relevant for 
both the classical multi-baker and the ID random walk, which is modeled by 
the classical system. Fixing the Planck constant (N = const) and performing 
long time limit we observe {(Ar) 2 (t)) = kt 2 /N, which is reflection of the 
crystal-like structure of the system. 

4-3 Extremal properties of the MSD 

To get another perspective on the above results consider a more general class of 
quantum multiplexer maps [17] where the local dynamics is given by a different 
map B, e.g a cat map, standard map, or other, replacing the baker map. Then 
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the formula (27) is still valid, however the eigenvalues and eigenvectors are 
those of the new map B. 

Using formula (27) it is easy to see that the time-dependent mean square 
displacement for any quantum multiplexer map has to satisfy 

(1) <(Arf(0)> = 

(2) <(Arf(l)> = 1 

(3) < <(Ar) 2 (t)> < 

In fact, those results are true in both the quantum and the classical case. More- 
over, it is possible to find local dynamics which realize both of the extremal 
cases. 

To see it, take as the local map a right-left exchange operator, defined classi- 
cally by 



B(n,x,y) :-- 



(n,x + l/2,y), for0<x<l/2, 
(n,x-l/2,y), for 1/2 < a; < 1. 



and quantum mechanically by 



B 



1 

1 



(using the context of Section 2). The dynamics is obvious: the particle jumps 
between two neighboring sites, leading to the mean square displacement having 
values for even and 1 for odd times. 

On the other hand, taking identity as the local dynamics, we induce ballistic 
motion: particle starting in the state going initially to the right will keep on 
going to the right, leading to the ballistic transport: 

<(Ar) 2 (t))=t 2 

Therefore, we see that translational invariance of the coupling operator T 
allows, in principle, for large asymptotic freedom: 

<(Ar) 2 (t)> oc f, 

where < a < 2 5 . An interesting question is, what can be realized in 



A more precise statement is: there exists a constant a £ [0, 2] such that 
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practice. While no constraints seem to be imposed on the classical level, the 
structure of the quantum mean square displacement, Eq. (26) suggests that 
for fixed H only a = or a = 2 are viable. One of the questions it raises is 
what conditions need to be imposed on the internal dynamics so that the 
semi-classical limit leads to anomalous diffusion a ^ 1. We interpret our 
RMT result as follows: internal fully chaotic dynamics (mixing) of the classical 
map implies generically diffusion. We are thus led to believe, that anomalous 
diffusion can arise in semi-classical limit in systems with partially chaotic, 
partially integrable internal dynamics. Similar observations were made often 
before in the context of different types of systems and transport in phase space 
as opposed to the transport in real space, which we discuss here. More precise 
results require further study. 



5 Comparison of numerical results with random matrix theory pre- 
dictions 

We have studied numerically a number of cases of quantum multi-bakers as 
well as some more general systems [17]. Here we present some of those results 
to show the quality of prediction afforded by the random matrix theory. 

Our numerical evaluation of the formula (27) shows that generically, almost 
every choice of phases defining the quantization (3) leads to results between 
the COE and CUE results. Figure 3 shows results obtained for the quantum 
multi-baker (2) with phases ip q = 0.73, ip p = 0.11. We show the mean square 
displacement as a function of time for N = 200 and for N = 2000 on log-log 
plot. There is a clear crossover between the diffusive and ballistic behavior 
on both plots. Figure (3.c) compares the values of J2j(\Jjj\ 2 ) calculated for 
the same quantum multi-baker for various values of N with the predictions of 
random matrix theory (45). We see that for iV > 50 the results lie between 
the analytic values for COE and CUE ensembles. 

In cases when the two phases add up to one tp q + tp p — 1, we find from our 
numerical work that the ballistic coefficient J2j{\Jjj\ 2 ) vanishes. Therefore, 
only the oscillating contribution remains in the expansion (26). We still ob- 
serve diffusive behavior up to Heisenberg time, although not as clearly as 
for generic systems. After this transient the particle localizes and the mean 
square displacement oscillates irregularly around some average value. Figure 4 
shows the mean square displacement for quantum multi-baker with phases 
cp q = cp p = 1/2 [9] for iV = 50 (top row) and N = 200 (bottom row). In 
both cases the left plot shows the short time behavior on the log-log scale. 
We see the same diffusive behavior as in the generic case. The right plots are 

lim^ 00 ((Ar) 2 (t))/t' 3 = oo for (3 < a, and lim^ 00 ((Ar) 2 (t))/t /3 = for a < [3 
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Fig. 3. Log-log plot of the mean square displacement in the quantum multi-baker (2) 
with phases ip q = 0.73, tp p = 0.11 as a function of time for N = 200 (top left) and 
for N = 2000 (top right). The crossover between the diffusive and ballistic behavior 
is clearly visible on both plots at the Heisenberg time t = N. The bottom figure 
shows the values of jj ^2j(\Jjj\ 2 } calculated for the same quantum multi-baker for 
various values of N with the predictions of random matrix theory (45). We see that 
for N > 50 the results lie between the analytic values for COE and CUE ensembles. 



normal scale and show the oscillations encountered on the long scale. All the 
plots are numerical evaluations of the exact formula (26). We suspect that 
this non-generic, localized behavior will be seen whenever the phases ip q , ip p 
sum to unity, but both the proof and a physical understanding remain to be 
developed. 



6 Multi-baker maps and quantum random walks 



We conclude this discussion of the quantum multi-baker maps by establish- 
ing connections between our work and some other research on quantum ran- 
dom walks. We restrict ourselves to discrete quantum walks on line. Con- 
tinuous quantum random walks [26,27,28] and discrete quantum walks on 
graphs [29,30] are less connected to our work. 

A popular model of a quantum random walk, known today as the Hadamard 
walk, was described some ten years ago by S. Godoy and S. Fujita [31]. It was 
obtained through approximating the evolution of some special wave packets in 
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Fig. 4. Mean square displacement for quantum multi-baker with phases 
if q = if p = 1/2 for N = 50 (top row) and iV = 200 (bottom row). The left plots 
show the short time behavior on the log-log scale. We see the same diffusive behav- 
ior as in the generic case. The right plots are normal scale and show the oscillations 
encountered on the long scale. All the plots are numerical evaluations of the exact 
formula (26). 

Kronig- Penney type potential. A similar model was discussed by D. Aharonov, 
et al. [32], who considered the motion of spin-| particles in one dimension, and 
proposed experimental realization of the walk in the framework of quantum 
optics. Further, D. Meyer studied quantum cellular automata for their possible 
applications to quantum computing. In section 5 of [33], Meyer obtained a 
quantum random walk model closely related to the Hadamard walk, discussed 
later by A. Ambainis et al. [16]. 

Looking for an analytically tractable model of quantum transport, we quan- 
tized [12] the classical multi-baker map, a deterministic model of classical ran- 
dom walk, obtaining a family of quantum random walk models, the simplest 
case being the Hadamard walk. We studied both the translationally invariant 
and random versions of the system. 

Since then a number of other papers appeared on the topics, e.g. [34,35,36,37,38,39,40,41,42,43,44,45]. 
Of particular interest for our work is the work of Brun et al. [46,47,48] who 
also studied the transition to classical behavior, but instead of considering 
semi-classical limit they assume interaction with environment implemented in 
decoherent coins. 

Much of the work on discrete quantum random walks can easily be refor- 
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mulated as a special case of the multi-baker maps (or multiplexer maps, see 
ref. [17]) in particular, the map with the largest possible value of Planck's 
constant, h = 1/2, or equivalently, N = 2. In order to make this connection 
we review, very briefly, the structure of quantum random walks as described 
by much of the literature. 



The simplest quantum random walk described in the literature is the Hadamard 
walk. Here one considers a collection of quantum states, which we will denote 
as \n,d). The parameter n will take on all integer values within a prescribed 
range, and here we take the range to be Z, all integers, positive and negative. 
The parameter d — r,l denotes a direction, to the right or to the left, respec- 
tively. Next, one defines an elementary quantum transformation, H, which 
takes | n, d) into symmetric and anti-symmetric combinations, as 



f \n,r')^ 
V K0 j 





(64) 



This step is now coupled to a translation, T of \n,r') one unit to the right 
and, similarly, \n, V) one unit to the left, as 



T\n, r) — \n + 1, r) 

T\n,l) = \n-l,l) (65) 

H is called the Hadamard gate in quantum computing literature [49,50], which 
gives rise to the name Hadamard walk. Finally the two operations are com- 
bined to form the operator T o H = W, which describes one step of the 
Hadamard walk, 

l,r,t-l) + |ra + l,Z,f-l)] 

l,r,t-l) -\n + l,l,t-l)\. (66) 

This is identical to Eq. (2) for N = 2 for the special case that the Balazs-Voros 
phases [8], tp q = ip p = 0, are used. Therefore, the multi-baker map and the 
Hadamard walk are identical, for the Balazs-Voros phases and N = 2. The 
Hadamard equations have been solved in detail for various boundary condi- 
tions [12,16,34,38,45]. Similarly, the N = 2 multi-baker map has been solved 
for a variety of phases and boundary conditions [12]. It should be appreci- 
ated that the multi-baker maps represent generalizations of quantum random 
walks in a number of respects. The inclusion of phases in the transformation 
equations allows one to treat a variety of maps including uniform, periodic, 
quasi-periodic, and random systems, something typically not considered in the 



\n,r,t) = — [\n- 
\n,l,t) = -y^[\n- 
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theory of quantum random walks. Thus one can find a range of phenomena 
in quantum multi-baker ranging from localization to ballistic motion. Such 
phenomena show up in many condensed matter systems. Further, the use of 
a variable Planck constant allows for several channels of motion to be taking 
place at once, and allows us to treat semi-classical as well as strong quantum 
versions of these maps. This is what makes the multi-baker maps so appealing 
for studying exact transport properties of condensed quantum systems. 

One common feature of different quantum random walks and quantum multi- 
baker maps is the ability to include various models of decoherence processes. 
In the multi-baker maps one should distinguish between quenched disorder 
and annealed disorder. In quenched disorder the phases tp q , tp p vary randomly 
for cell to cell, but once values are given for each cell, the phases keep those 
values for all time steps. This allows some memory to be generated by the 
dynamics and results in localization [12]. For annealed disorder, the phases 
are chosen randomly for each cell and at each time step. The phases also may 
be specified up to some additive noise contribution. These processes represent 
the effects of random external noise and result in some degree of decoher- 
ence [36,46]. Further, in both quantum random walk models and multi-baker 
maps one can produce decoherence by replacing the unitary matrix describ- 
ing the time evolution by a completely positive super-operator. This allows for 
both dissipation and decoherence. As this topic requires a separate discussion, 
we will not pursue it here, but leave the discussion of decoherence for further 
publications. 



7 Discussion of the results and other comments 

This paper described calculations of important properties of the quantum 
multi-baker map, a simple model for transport in quantum systems. This map 
is the quantum version of a well known, classically chaotic system. In partic- 
ular we were able to show, using random matrix theory, that it is possible to 
obtain an analytic expression for the equilibrium mean square displacement 
of a particle whose motion is governed by this quantum map. This expression 
is of particular importance for an understanding of transport in quantum sys- 
tems because it allows us to carefully examine the particle's average motion 
both for finite times and for non-zero values of Planck's constant. Thus we can 
see how the two non-commuting limits, t — > oo or h — > 0, interact with each 
other. While our results are not surprising, in the sense that we already knew 
qualitatively the limiting forms, we did not have a general expression for the 
mean square displacement even for simple quantum systems before this work 
was carried out. Moreover, the results from random matrix theory agreed very 
well with careful numerical evaluations of the mean square displacement for 
the multi-baker map. This agreement is also very important because random 



27 



matrix theory does not really apply to any individual system, as far as one 
knows, but rather to the average behavior of an ensemble of similar systems. 
Thus in the absence of a rigorous justification of the use of random matrix 
methods, their predictions must always be checked by comparisons with nu- 
merical studies. 



It will be interesting to find other signatures of classical behavior in the quan- 
tum multi-baker map that are more closely connected to the chaotic behavior 
of the classical map than the mean square displacement. For example, is there 
any trace of the classical, hyperbolic dynamical behavior in the quantum ver- 
sion of the map? One prominent feature of the classical multi-baker model 
is the formation of fractal structures on arbitrarily fine scales as an initial 
non-equilibrium distribution of points relaxes to a final uniform equilibrium. 
This fractal structure can, in turn, be related to a positive entropy production 
typical of the approach to equilibrium of a macroscopic system. The quan- 
tum version of this phenomenon is largely unexplored and is of considerable 
interest to us. 



Further work will also be devoted to exploring the consequences of the free- 
dom one has in choosing phases of the unitary operators. This freedom allows 
us to consider random, quasi-random, and periodic systems, each with its 
own interesting physical properties. The quenched random system has been 
explored to some extent [12], and provides a simple model of localization in 
one-dimensional disordered systems. Here the classical limit is expected to be 
recovered through a phenomenon where the localization length grows infinitely 
large as Planck's constant approaches zero, such that almost normal diffusion 
takes place within each localized region. This picture needs to be verified, and 
further systems still need to be explored. 



We also began to establish connections between this work on quantum multi- 
baker maps and a body of related work on quantum random walk processes. 
Much of that literature is, in fact, devoted to the simplest multi-baker model, 
although this connection has not been recognized in the random walk litera- 
ture. Here we showed that the N = 2 quantum multi-baker process is identical 
to the Hadamard walk of random walk theory. Although quantum computing 
is still far from being realized [49,50], it is possible that the theory of quan- 
tum random walks may be of some value for constructing algorithms based 
on quantum walks that will be faster than current classical random walk al- 
gorithms [27,43]. Future work will be devoted to expanding the connections 
between the quantum random walks and multi-baker processes for a variety 
of model systems. 
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